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We investigate a superconducting quantum interference device (SQUID) based on carbon nan- 
otubes in a fork geometry [J. -P. Cleuziou et ai, Nature Nanotechnology 1, 53 (2006)], involving 
tunneling of evanescent quasiparticles through a superconductor over a distance comparable to the 
■ superconducting coherence length, with therefore "non local" processes generalizing non local An- 

»' ' dreev reflection and elastic cotunneling. Non local processes induce a reduction of the critical current 

, and modify the current-phase relation. We discuss arbitrary interface transparencies. Such devices 

in fork geometries are candidates for probing the phase coherence of crossed Andreev reflection. 
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I. INTRODUCTION 



Implementing experimentally!^^ and understanding theoreticallj^^>Sii8i9>ioai>i2aaia5a6^a8a9.20^.22^i24^ the 
' possibility of emitting spatially separated electron pairs from a superconductor in different electrodes has aroused 
, considerable interest recently, in connection with the realization of a source of entangled pairs 

of electrons21iS2. The 

^ ' Josephson junction through a carbon nanotube quantum dot^^ and the superconducting quantum interference device^l 
QJ (SQUID) realized recently can be viewed as steps towards future implementations of transport of spatially separated 
pairs of electrons in quantum information devices based on carbon nanotubes^^. As another application, the carbon 
nanotube SQUID realized recently by Cleuziou et al^ in the fork geometry in Fig. [T] proves the feasibility of future 
. measurements of magnetization reversal of individual molecular magnets. 

Even more miniaturized devices may be realized in the near future, with geometrical dimensions comparable to 
an intrinsic length scale of the superconductor: the characteristic length?^ associated to the superconducting gap 
1^0 1- Such devices^^i^liSS may be sensitive^i^ to the fact that Andreev reflection^ takes place in a coherence volume 
C of linear dimension therefore allowing for the possibility of splitting Cooper pairs in two parts of the circuit. 
^ ' Andreev reflection^ is the process by which a spin-up electron incoming from the normal side on a NS interface 
between a normal metal N and a superconductor S is reflected as a hole in the spin-down band while a pair of 
electrons is transmitted in the superconductor. In a NaSN;, structure with the electrical circuit on Fig. [21 an electron 
in electrode Na can be scattered as a hole in Nf, if the contacts are separated by a distance comparable to the 
superconducting coherence length (see Fig. [^h- for non local Andreev reflection). Alternatively, an electron from 
Na can be transmitted as an electron in N;, across the superconductor® (see Fig.l^b for elastic cotunneling). 
\ The goal of our article is to address possible realizations of non local Andreev reflection in future carbon nanotube 
' SQUIDs. These devices^i^^ would provide further— experimental signatures of the phase coherence of Cooper pair 
splitting. Compared to the previous Refs. [24.25] we investigate here higher order processes in the tunnel amplitudes 
giving rise to "non local Andreev bound states" . By contrast, if the distance between the Josephson junctions is much 
— . _ larger than the superconducting coherence length, the bound states are "local", and localized over two separated 
, regions of extend on each Josephson junction, not coupled by non local processes through the superconductor S 
L; ; (see Fig. HD. 

. ^ ' The article is organized as follows. Preliminaries are presented in Sec. |TT1 Our results are presented in Sec. IIII AI 
for the dc- Josephson effect in single channel systems. Multichannel effects are discussed in Sec. IIII Bl Concluding 
H remarks are presented in Sec. IIVI Some details are left for Appendices. 

_ _ I 

II. PRELIMINARIES 



A. Carbon nanotube superconducting interference device 



A superconducting quantum interference device (SQUID) made of two superconductors and a carbon nanotube 
is considered (see Fig. [1^), according to the recent experiment by Cleuziou et al^. As a natural hypothesis, the 
proximity effect between the superconductor and the carbon nanotube (i.e. the penetration of pairs from the super- 
conductor to the nanotube) is supposed to induce a minigap |Ao| in the portions of the nanotube in contact with the 
superconductors. 
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The nanotube is divided in five sections connected to each other, from top to bottom: superconducting top section 
with a minigap |Ao| in contact with the superconductor S'; quantum dot number 1; superconducting middle section 
of the nanotube with a minigap |Ao| in contact with S; quantum dot number 2; and superconducting bottom section 
with a minigap |Ao| in contact with the superconductor S' (see Fig. [T^). 
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FIG. 1: (Color online.) Schematic representation of the carbon nanotube SQUID (a). The hopping description in the ofT- 
resonant state is shown on (b). The transport properties of the SQUID depend only on the enclosed flux and on the su- 
perconducting phase differences, (c) shows a double bridge between two superconductors, (d) is a double insulating bridge 
between two superconductors, made of two multichannel insulators in parallel. Only non local processes through one of the 
superconductors (the superconductor S) are allowed in all these situations. The Oxyz axis is shown on (b). 



Depending on the value of the gate voltage in experiments, the quantum dots 1 and 2 can be tuned from off-resonant 
to resonant (changing the gate voltages has the effect of shifting the dot energy levels). Our article discusses mostly 
the off-resonant state (in short: off-state) such that dots 1 and 2 have a vanishingly small density of states within 
the minigap (see Fig.[3K). It was well established experimentally by Cleuziou et alM- that their SQUID can be tuned 
from the off-state with a very small critical current to the on-state with a large critical current by changing the gate 
voltages coupled to the two quantum dots formed by portions of the nanotube in between a and a', and in between 
(3 and fi' (see Fig.[I]a). 

Our modeling is intended to capture "non local bound states" involving multiple electron-hole processes back 
and forth between dots 1 and 2. We make the following simplifying assumptions. First, proximity effect between 
the nanotube and the superconductor is not described explicitely on a microscopic basis: for highly transparent 
interfaces between the carbon nanotube and the superconductor S, we treat proximity- induced superconductivity in 
the nanotube as bulk superconductivity, and therefore we use the denomination "gap" instead of "minigap" . Non 
local processes then take place at the discontinuities of the superconducting order parameter at a and a' . For lower 
interface transparencies between the superconductor S and the carbon nanotube, electrons and holes may propagate 
in the portion of the nanotube in contact with S before undergoing non local Andreev reflection or elastic cotunneling, 
which amounts to averaging over many channels for non local processes. 

Second, we consider the off-state with a vanishingly small density of states if the absolute value of energy is 
smaller than the gap (see Fig. [3^) , and we use a standard description as tunnel amplitudes^i^^i^ connecting the 
superconductors S and S'. 
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(a) Crossed Andreev reflection 




(b) Elastic cotunneling 



FIG. 2: (Color online.) Schematic representation of the electrical circuit for probing the non local conductance in normal metal 
- superconductor - normal metal (NaSNj,) structures'^, (a) shows a schematic representation of the non local Andreev reflection 
process changing a spin-up electron in electrode N;, into a spin-down hole in electrode Na and leaving a Cooper pair in the 
superconductor, (b) is a schematic representation of elastic cotunneling transferring an electron from one normal electrode to 
the other via a trip through the superconductor. 
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FIG. 3: (Color online.) Schematic representation of quantum dot density of states, with a dot level spacing 5 and a level 
broadening F. (a) corresponds to the off-resonant state considered in our article, with no resonant level in the gap window, (b) 
corresponds to the resonant situation with a large density of states within the gap. 
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As a third assumption, it is well known that a single wall carbon nanotube contains two conduction channels. The 
idealized case of a single hopping amplitude is discussed in Sec. IIII Al Multichannel contacts are left for Sec. IIII Bl 

As a final assumption, Coulomb interactions in the quantum dots are not accounted for, so that S on Fig. [3] is 
supposed to be a finite size effect not due to Coulomb interactions. 

The distance Ra,f3 between a and f3 (see a and /3 on Fig. [T]d) is supposed to be comparable to the coherence 
length ^0- The shortest path connecting a' and /3' (see a' and /?' on Fig. [1)3) is much larger than {Ra',0' ^ (,q), 
with therefore no "non local" quasiparticle tunneling between a' and (3' . Compared to the fork in the experimental 
geometry realized in Ref. (27| , these assumptions on geometry can be implemented in future experiments by reducing 
the distance between a and jS (see a and /? on Fig. [T)d) down to values comparable to the superconducting coherence 
length ^0- In the experiment by Cleuziou et al^, the distance i?a,/3 between a and f3 is comparable to the distance 
Ra,a' between a and a', and to the distance -R/3,/3' between /3 and /3' (of order 400 nm). We consider on the contrary 
future fork geometries with Ra^/j ^ and with Ra^p ^ Ra.a', Rp,i3'- 

Choosing the gauge = —By/2, Ay = Bx/2 and A^ = Q (with the Oxyz axis on Fig. [1)3), with B the apphed 

magnetic field and A the vector potential, leads to a finite value for Adr, and to Adr = — Jp Adr if we 
suppose Ra.a' = -R/3,/3'- Givcu Ra,i3 ^ Ra,a' , wc ncglcct the line integral of the vector potential between a and f3: 
Adr ~ 0. The line integral of the vector potential along a path from a' to f3' in S' is finite but quasiparticle 
propagation from a' to p' has a vanishingly small amplitude (because Ra'^p' S> ^o)- Non local processes between a' 
and /?' are thus negligible. 



The SQUID corresponds to two Josephson junctions, one for each interface. If the junctions are far apart (at a 
distance much larger than the superconducting coherence length), and for single channel weak links, one negative and 
one positive energy Andreev bound state is located on each junction, therefore leading to a total of four Andreev 
bound states for the SQUID (two bound states at positive energy with respect to the Fermi level, and two bound 
states at negative energy). As we show below, non local processes induce a coupling between these bound states in 
the form of level repulsion. 

The supercurrent is obtained from differentiating the free energy with respect to the superconducting phase differ- 
ence Aip. Beenakker'^*' finds three terms contributing to the supercurrent, some of which date back to the early stages 
of Josephson junction theory'^^. First at zero temperature the following term corresponds to a summation over the 
discrete bound states within the gap: 



where (p is the flux enclosed in the loop of the SQUID, and where Nabs is the number of Andreev bound states. 
The step function in energy 9 [— J7„(A(/j, (j>)] selects Andreev bound states below the Fermi level. The second term in 
Ref. [36], corresponds to the contribution of the continuum to the supercurrent. The contribution of the continuum is 




not included in the discussion in the forthcoming Sees. IIII Al and ITlI Bl We will justify in Sec. IIII A^ that it is indeed 
negligibly small in the situations that we consider. The third and last term in the expression of the supercurrent 
obtained by Beenakker— vanishes if the superconducting gap is independent on the phase difference, which we assume 
in the following. 

The bound states are obtained in a standard description as the poles of the fully dressed Green's functions. The 
later is determined by the Dyson equations, which allows to describe weak links ranging from tunnel contacts to highly 
transparent interfaces. 

The Green's functions of the superconductor are obtained by Fourier transform in a well known procedure^. For 
superconductors isolated from each other, the local Green's function takes the form 



B. Microscopic Green's functions 




(1) 
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where the superconducting phase variables Lpi^ and Lpji — (fL + are shown on Fig.[T|3, jAoj is the superconducting 
gap and pjy the normal density of states. 
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The non local Green's functions take the form 

5a,/3H = 90A^)=Cic.)7rp^l f IA.I.;!L-,..A )cos(fc^^i?„.,) 



^ 1 Sin(fci;^i?a,^) i- , (4) 



where we use the notation C(a;) for 

C{lo) ~ exp 



(5) 



We parameterize below the strength of non local processes by Cq = C{uj = 0) = exp [—2Ra^p/^o], with = = 0). 
The strength of non local processes is parameterized by Cq, ranging from an absence of non local processes (Cq — 
for Ra,i3 ^ ^o) to non local processes taking their maximal value (Co — 1 for Ra^f3 <C ^o)- One has then the following: 

C{lu) = [Caf^^^^° . (6) 

Intermediate values of Co are expected for a carbon nanotube with Ra^is of order ^o- Note that in the case of 
three dimensions (not considered here), the cos (kpRa^p) and sm{kpRa^f3) factors are interchanged, and C{uj) — 

exp [-2Ra,0/^{Lj)]/Ra,i3. 

The condition Ra',p' » Co (see Fig. [TJd) leads to 

9a',0'{(^) = 9/3', a' (OJ) =(^00^- (7) 

The fully dressed Green's functions at energy w are obtained via the Dyson equations taking the following form in 
a compact notation: 

G{uj)^g{uj)+g{uj)(E}tt(E>G{uj), (8) 

where g(w) corresponds to the Green's functions of the superconducting electrodes isolated from each other in the 
absence of tunnel amplitudes, St is the Nambu hopping self-energy, and G{oj) is the fully dressed Green's function. 
The notation ® denotes a convolution over the network labels a, /3, a' and /3' (see the notations on Fig. [TJa). For 
instance one has the following: 

Ga,l3iuj) = ga,p{uj) + ga,a{uj)ia,a'Ga',f3i(^) + ffa,/? (a;)f,3,;3' G^g' ,/3 (w) ■ (9) 

The set of fully dressed Green's functions arc then obtained from matrix inversion, and the Andreev bound states 
correspond to the poles within the gap. They are determined either from the corresponding analytical expressions of 
the Green's functions, or from a numerical solution. 

III. RESULTS 
A. DC Josephson effect for a single transmission channel 

1. Amplitude and minima of the critical current 

We find a reduction of the supercurrent upon increasing the strength Co of non local processes (see Fig. |4|) . On 
this figure, the critical current is averaged over all realizations of the Fermi phase factor kpRajj corresponding to 
averaging over many samples with different Fermi phase factors. The interfaces of a superconducting electrode are 
not controlled on atomic scale and it is thus a natural assumption^'^"* to use a uniform distribution of the Fermi phase 
factors kpRa^0. As expected, the reduction of the supercurrent by non local processes in a collection of single channel 
systems is in agreement with a collection of multichannel systems (see Sec. IIII B|) . 

Anticipating the forthcoming Sec. IIII A^ we note that non local Andreev reflection changes an electron with positive 
energy on one junction into a hole with negative energy on the other junction. As a consequence, bound states with 
positive energy are coupled to bound states with negative energies. The resulting level repulsion among bound states 
with opposite energies reduces in absolute value the slope of the phase dependence of the Andreev levels and therefore 
reduces the critical current. We evaluate in Appendix IB] the critical current for tunnel interfaces, and we confirm by 
this analytical treatment the reduction of the critical current by non local processes. 
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FIG. 4: (Color online.) Critical current as a function of the magnetic flux in the loop of the SQUID, with {Ta,Tp) — (0.8, 0.8) 
(a) and with {Ta,Tp) = (0.8, 1) (b). The critical current is averaged over all values of kpRa,!!- Different curves correspond to 
different values of Co, the strength of non local processes. 



2. Level repulsion among Andreev bound states and current-phase relation 



Assuming a distance between the Josephson junctions much larger than the coherence length, and assuming also 
single channel contacts, we find one Andreev bound state with negative energy localized on each junction, as it 
should. The bound states extend in the superconductor over a region of size comparable to the coherence length. If 
the distance between the Josephson junctions becomes comparable to the coherence length, the bound state energy 
levels depend on the coupling corresponding to "non local" propagation in the superconductors (see Fig.[T]). 

The variations of the bound state levels (±i7i, ±1^2) with the flux cj) enclosed in the loop are shown on Fig. [5^ for 
Co = (absence of non local processes) and on Fig. [SJs for Cq = 1 (maximal value of non local processes). 

Level repulsion among Andreev bound states upon increasing Co (see Fig. [5}d) has the effect of reducing the slope 
of the bound state energy levels versus phase relation, which reduces the supercurrent. The bound states take the 
simple form given in Appendix [Cl (see Eqs. (jCip and (|C2[) ) for a symmetric contact with ta = tt- We do not present 
in the article the too heavy expression of the bound state levels in the general case. 

Now we consider single channel transmission modes between the superconductors S and S', and with non local 
processes at the interfaces where a step function variation of the superconducting gap is assumed. As seen from 
Fig. [6l the SQUID current phase relation fluctuates from sample to sample. 

We conclude this section by discussing the contribution of the continuuni^^ for the hopping model of SQUID in the 
off-state. The supercurrent is obtained as the integral over energy of the spectral supercurrent. We show on Fig. [7] a 
typical variation of the spectral supercurrent as a function of energy. We find practically no contribution of energies 
larger than |Ao| in absolute value, as opposed to other cases such as Ref. [3!||. We conclude that the contribution 
of the continuum is negligible for the hopping model of SQUID in the off-state in which the hopping elements are 
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FIG. 5: (Color online.) Variation of the bound state levels in the absence of non local processes (Co = 0)(a) and with the 
maximal strength of non local processes (Co = 1) (b). The bound states repel each other for Co = 1 in the presence of non 
local processes on (b). For this figure, we use Aip = 0, kpRa.p = 1 + 2nn (with n an integer), Tq, — 0.8 and Tf; = 1. A similar 
repulsion between Andreev bound states is obtained in the dependence of the bound state energy levels as a function of the 
phase difference A 95. 



energy-independent. The supercurrent is therefore well approximated by Eq. ([T]) as in Sees. IIII Al and IIIIBI As 
a physical interpretation, Andreev bound states are localized in the superconductor in a region of size set by the 
coherence length. A single channel weak link coupling two superconductors^ is a very localized perturbation which 
does not couple most of the extended states in the superconducting electrodes, which explains why the states of the 
continuum are almost insensitive to the phase difference between the superconductors in the considered geometry 
with localized interfaces. 



B. SQUIDs involving multichannel contacts 



A metallic carbon nanotubc consists of two conduction channels, and it is thus natural to extend the discussion in 
Sec. IIII Al to the case of multichannel^'^ junctions (see Figs. [TJ: and d). We evaluate the density of Andreev bound 
states^'"' for multichannel systems with Nch = 2 and Nch = 15 channels, and with Cq = and Cq — 1 (see Fig. [8|). 
Increasing the strength of non local processes by increasing Cq reduces the density of Andreev bound states, therefore 
reducing the value of the supercurrent, in agreement with Sec. IIII A II 



IV. CONCLUSIONS 



To conclude, we have discussed signatures of non local Andreev reflection on the current-phase relation of a dc- 
SQUID in a fork geometry similar to Cleuziou et alM; but with the dimension of the middle superconductor compa- 
rable to the superconducting coherence length, so that quasiparticles may tunnel through the superconductor, with or 
without electron-hole conversion. Compared to a geometry consisting of two parallel normal bridges connecting two 
superconductors'"*'^^, we investigated here processes of higher order that are not washed out by disorder. For idealized 
single channel systems with sharp step-function variations of the superconducting gap in the nanotube, we found that 



8 



Supercurrent vs phase difference 



CM 



< 
CM 



CM 




kpR= 
kpR= 



0.0 -H — 

:0.5 -X — 
:1.0 >K 



-0.25 0.25 
A(p/27t 

kpR=1.5 B 



kpR=2.5 



kpR=3.0 A 
k^R=3.5 



FIG. 6: (Color online.) Supercurrent versus phase difference with non local processes, for different values of kpRa^p, and 
for {TcTjj) = (0.8,0.8), without magnetic flux (j> = (a), with magnetic flux <j) — 2 (b). The supercurrent is very small for 
kpRa,/! ~ 1.5 + 27rn (with n and integer), not far from kpRa^/s — n/2 + 2nn. 



non local processes induce sample to sample fluctuations of the current phase relation due to the dependence of non 
local processes on the Fermi phase factors. Multichannel systems capture moderate interface transparencies between 
the nanotube and the superconductor because in this case electrons incoming in the superconductor can propagate 
in the nanotube before undergoing crossed Andreev reflection. Increasing the strength of non local processes reduces 
the supercurrent, as for single channel systems. 

From the point of view of future experiments, non local processes play a role in SQUIDs with fork geometries and 
with junctions made of carbon nanotubes, normal metals or semiconducting quantum wires. It would be interesting 
to measure the reduction of the current-phase relation of the SQUID upon increasing the strength of crossed processes 
in multichannel systems, via a comparison of samples with different dimensions, or via the temperature dependence 
of the coherence length. A more difficult experiment consists in probing sample to sample fluctuations of the SQUID 
supercurrent. A good characterization of the sample parameters (such as number of channels, interface transparencies) 
is required in order to distinguish between the intrinsic fluctuations of crossed processes and unwanted variations of 
the junction parameters when changing from one sample to another. As pointed out to us by F. Giazotto, the 
strength of non local processes can be monitored by the temperature dependence of the superconducting coherence 
length ^0 = fivp/l^ol (in the ballistic limit) or ^0 = -y/ HD / \ Aq | (with D the diffusion coefficient in the diffusive 
limit) because the superconducting gap decreases with increasing temperature. Increasing temperature has thus the 
effect of reducing |Ao| and enhancing the coherence length. It is expected that the total amplitude of supercurrent 
decreases with increasing temperature, but the relative contribution of non local processes increases upon increasing 
temperature. 
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FIG. 7; (Color online.) Variations of the spectral supercurrent for a given realization of the microscopic Fermi phase factors in 
the vicinity of hui = huj/2'K ~ — |Ao|. We find practically no contribution of the continuum for the hopping model of SQUID 
in the off-state. We use the phase difference = 2 on the figure, but similar variations of the spectral current are obtained 
for other values of Aip. We also chose Co = 1 on the figure. 



APPENDIX A: FLUCTUATIONS OF NON LOCAL TRANSPORT THROUGH A DIFFUSIVE 

SUPERCONDUCTOR 

The charge transmission coefficient of a diffusive superconductor vanishes after averaging over disorder in the 
diffusive limitiiii^. To show that the charge transmission coefficient of a disordered superconductor fluctuates at 
the scale of the Fermi wave-length Xp, we show that the opposite hypothesis does not hold. Simply, we obtain 
fluctuations of the charge transmission coefficient by adding a very small extra ballistic region much larger than the 
Fermi wave-length but much smaller than the elastic mean free path. 



APPENDIX B: CRITICAL CURRENT IN THE TUNNEL LIMIT 



In this Appendix, the supercurrent is expanded in the tunnel amplitude ta.b connecting the two interfaces between 
the superconductors. We do not detail the corresponding calculation based on diagrammatic perturbation theory. We 
start with the first terms of an expansion of the supercurrent in the tunnel amplitudes [see the following Eqs. (jBip - 
(jBlOp j. The dimensionless parameters Ta.b — {T^PNtafiY , related in the tunnel limit to the dimensionless interface 
transparencies through Ta^ — ^Ta.bi are supposed to be small. The Dyson equations are expanded systematically 
in the tunnel amplitudes and the lowest order diagrams are collected, leading to the following expansion for the 
supercurrent: 

/s(Av3, 0) = sin (A(/5 + 0) + A sin (Av3 - (/)) -I- B sin (2 A(/3) (Bl) 
-I- i:'aSin(2(A(/3 + 0))+DbSin(2(A(/7-0)), 

with 

Aa = -^T.IAol (B2) 
h 

Ab = -^^blAol (B3) 
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FIG. 8: (Color online.) Distribution of the Andreev bound state density of states at phase difference Ay = vr, for Nch ~ 2 
channels (a) and for Nch ~ 15 channels (b), without non local processes (Co ~ 0, □) and with non local processes to their 
maximal value (Co = 1, ■). Each histogram is normalized to the number of Andreev bound states averaged over the realizations 
of the Fermi phase factors. 



B = -—TauC^JAolsin^ {kFRc.,p) (B4) 

Da = -|^^a|Ao| (B5) 

D, - -|^^b|Ao|, (B6) 

where C, is the value of the parameter C{uj) where huj takes the value of the bound state energy level. In this case, 
the bound states are very close to the gap, so that C* ~ 1. The contribution of non local Andreev reflection to the 
supercurrent (term B in Eq. (|B4[) ) is, in the tunnel limit, much smaller than the contribution of local tunneling of 
Cooper pairs from one superconductor to the other (terms Aa and At in Eqs. (|B2|) and (|B3|) ). 

The reduction of the supercurrent upon including non local processes is described in the tunnel limit by an expansion 
of Eq. (|Bip around — Tr/2: (j) = 7r/2 + 5(j). We consider an ensemble of single channel systems with a collection of 
kpRa./s and thus we average to 1/2 the factor sin^ {kpRa.p) in Eq. (|B4[) . We define Ai^s^"^ as the value of /S.ip such 
that 

c)/5(Av. = Ay>(°),^) 

d{K^) = °' (^^^ 

and we obtain 

(0), (1 - t)54> ± V((l - t)5^? + 2r^(2 - 01/2^ 
cos ( A.," ) ^ 2r(2-Cg/2)^ ' ^^^^ 
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where we supposed a symmetric contact with Ta,b = For S4> ^ t, the critical current is given by 

IciTT/2 + Sc^)^iTTT\Ao\S(j), (B9) 

and for Scj) <^t, it is given by 

/,(7r/2 + 50) =^t2|Ao| (^2-^^ . (BIO) 
For = the expansion of the critical current is given by 

where non local effects do not enter Eq. (jBlip at the leading order. We conclude in the tunnel limit to a reduction of 
the contrast of the critical current oscillations upon increasing the strength C, of non local processes. The main body 
of the article corresponds to interfaces with moderate of large transparencies, described by Andreev bound states 
obtained from the Green's function dressed by tunnel processes to infinite order, as opposed to the limit of tunnel 
contacts considered in this Appendix. 

APPENDIX C: ANDREEV BOUND STATES 

In a symmetric SQUID with ta ~ h and no magnetic field (0 = 0), the bound states take the form 



nfiA^) = ±|Ao|v/^+(A^)/v/i?+(A^) (CI) 



nfiAifi) = ±|Ao|V^-(A^)/v/B-(A(^) (C2) 

with 

A^{Aip) = 2cos(A(^)T(l ±asin(/cFi?a./3)) + 1 (C3) 
(1 ± sm{kFR„^[^)f + T^Cl cos(/cF^a,/3)^ 

and 

B^{Aip) = [l + T{l±C.,sm{kpRa^fi))f +T^Clcos{kFRa,f3f (C4) 
where t = Tr^p^t^ = t'^ /W'^, related to the normal transmission by the relation^ 

nN^«JW')/i^ + tlb/WY (C5) 

where W = l/irp^ is the band-width, with pj^ the normal density of states. The bound state levels are determined 
self-consistently in such a way as C, is the value of C{uj) [see Eq. ([S])], where Hcu is replaced by the bound state energy 
in a self-consistent manner. 

Lets us consider the case kpR^^f^ = 27rn (with n an integer). The bound states levels deduced from Eqs (|C3p and 
(IC4[) are then degenerate: 



17J-(A^) = f72"(A(p) = -|Ao|Vl^^Q^shI(A^72)2 (C6) 

with 

» = 4./((l + + r^Cl) ^ ^,^,.|^l[7lcltVW' - 

where we introduced the band- width W according to Ref. [s^l ■ The SQUID is then equivalent to two identical S-I-S 
junctions, as seen from comparing Eq. (|C6p to Ref. (3^ . 

For kpRajj = 7r/2 + 27rn (with n an integer) the degeneracy is removed only if Co ^ 0: 



n^^i^^p) = -|Ao| Jl - a(i'2) sin(A(p/2)2, (C8) 
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with 



a 



(1,2) _ 



4t(1 ± C* 



(1 + t(1±C,))2 



(C9) 



where the "+" and "— " signs correspond to "1" and "2" respectively. 

In the presence of a magnetic field 7^ 0), the bound state energy levels take the form 



with 



f7j;_2(A^) = -|Ao|\/l - a(i^2) sin(A<y9/2)2 - /3(i,2) sin(A^) (CIO) 



(1,2) ^ 4TCOs(20)(l±a) 

" (l + T(l±C))2-4rsin(0)2(l^a)±4T2C,sin(20)2' ^ 

^(1,2)^ 2Tsin(2c/.)(lTa) ,^ . 

^ (l + r(l±C))2-4Tsin(0)2(l^C*)±4r2C*sin(20)2- ^ > 



and 



lA I lA I 1 ,_ sinW^(lTC)±rC.sin(20)2 

|Ao| = lAolWl - 4r + ±c,))2 " ^^^^^ 
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